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Abstract 

Generalized  parameters  characterizing  the  state  of  the  interplanetary  me¬ 
dium  (IM)  include  the  functional  form  and  the  rate  of  decline  of  charged  particle 
flux  in  solar  energetic  particle  (SEP)  events.  The  shape  of  the  particle  flux  decline 
is  of  particular  importance;  power-law  time  dependence  indicates  the  dominance 
of  diffusive  propagation,  whereas  exponential-law  decline  emphasizes  adiabatic 
deceleration  and  convection.  A  statistical  investigation  extended  for  a  long  period 
suggests  that  about  90%  of  SEP  decays  are  characterized  by  exponential  declines. 
Variation  of  the  decay  time  with  energy  and  angular  distance  of  the  observer  from 
the  flare  are  considered  as  well. 

1.  Introduction 

This  paper  considers  SEP  events  after  their  maxima  having  exponential 
decay  phases  between  1973  and  2001  as  a  continuation  of  a  systematic  statistical 
study  of  decay  phases  of  SEP  events  [1,2].  Here  we  analyze  the  dependence  of 
T  on  the  flare  site  and  on  proton  energy  for  all  events  satisfying  our  selection 
criteria.  Distributions  of  durations  of  shock-associated  and  non-associated  events 
are  also  considered. 

2.  Distributions  of  total  durations  of  decays  and  their  functional  form 

r- 

Time  intervals  were  selected  on  the  basis  of  a  homogeneous  data  set  ob¬ 
tained  by  the  CPME  instrument  aboard  IMP  8  between  1974  and  2001  [1,2]  during 
the  decay  phase  of  SEP  events  following  flares,  CMEs  and  interplanetary  shocks  of 
different  origin.  Only  those  parts  of  declines  were  examined,  which  could  reason¬ 
ably  be  described  by  exponential  dependence.  Fig.  1  displays  the  distributions 
of  the  total  durations,  AT,  of  shock-associated  exponential  declines  and  those 
without  shocks.  The  distributions  differ  insignificantly  for  AT  >  10  hrs,  whereas 
for  shorter  events  the  relative  number  of  shock  associated  events  is  30-40%  higher 
than  without  shocks.  This  can  be  due  to  the  fact  that  at  our  selection  criteria 
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short  duration  events  often  represent  regions  of  shock-trapped  particles,  which 
rapidly  pass  by  the  observer.  In  most  cases,  when  the  time  profile  following  the 
maximum  is  smooth  and  free  of  disturbances,  the  exponential  function  gives  the 
best  approximation  in  nearly  90%  of  the  cases  for  proton  energies  above  1  MeV. 


Fig.  1.  The  distribution  of  decay  times.  Fig.  2.  The  distribution  of  a  values. 

The  comparison  of  experimental  values  of  characteristic  decay  times,  r, 
with  those  obtained  previously  in  different  theoretical  models  [4,6]  shows  that 
theoretical  r  values  are  reasonably  close  to  the  fitted  slopes  in  nearly  50%  of  all 
cases  if  one  uses  the  average  F.solar  wind  speed  values  measured  later  when  the 
corresponding  plasma  in  which  particles  were  convected  arrives  to  the  observer. 
This  might  be  a  surprisingly  high  number  bearing  in  mind  that  in  rare  case 
remains  V  constant  for  a  sufficiently  long  time. 

3.  Energy  dependence  of  characteristic  decay  time  in  exponential  de¬ 
cays  of  proton  SEP  events 


To  study  the  energy  dependence  we  selected  events  with  exponential  de¬ 
cays  for  which  r  could  be  determined  within  an  accuracy  of  10%  for  high-energy 
differential  channels  (4.6-48  MeV),  i.e.  sufficiently  powerful  events,  and  where  the 
time  profiles  following  the  maxima  were  smooth  and  free  of  considerable  distur¬ 
bances.  Here  the  fluxes  of  4.6-15  MeV  proton  usually  reach  1  p/(cm^  s  sr  MeV). 
109  such  events  were  found  during  the  period  1973-2001,  and  the  exponent  as¬ 
suming  power-law  energy  dependence  r  =  CE-°‘  was  determined  (E  is  the  proton 
kinetic  energy).  Due  to  the  wide  energy  intervals,  the  errors  in  determining  r  was 
about  10%  due  to  not  strict  exponentiality  and  small  flux  variations.  The  accu¬ 
racy  of  determination  of  a  is  about  30%.  Presented  in  Fig.  2,  the  distribution  of  a 
values  can  be  considered  to  consist  of  three  separate  groups;  (1)  r  independent  of 
proton  energy  (-0.1  <  a  <  0.1,  40  events);  (2)  r  decreasing  with  energy  (a  >  0.1, 
57  events);  and  (3)  a  small  group  with  r  increasing  with  energy  (a  <  —0.1,  12 
events).  Thus  predominantly,  r  is  either  independent  of  or  decreases  with  proton 
energy.  A  constant  r  results  in  spatial  and  temporal  invariance  in  the  spectra  of 


energetic  paxticles  in  gradual  solar  events  discussed  in  [7]  when  nearly  identical 
spectra  are  seen  over  large  heliolongitude  region  some  time  before  or  after  shock 
arrival  depending  on  the  relative  positions  of  the  flare  and  observation  point.  We 
assume  that  it  is  possible  to  interpret  the  second  group  (decreasing  r)  as  having 
a  ’’memory”  of  diffusion  propagation  at  the  early  stage  near  the  Sun.  Indeed, 
diffusional  decay  time  profile  (power-law  shaped)  is  proportional  to  where 

D  =  Xv/ 3  (A  is  scattering  mean  free  path,  v  is  the  particle  speed)  is  diffusion  co¬ 
efficient  increasing  with  the  particle  energy  under  reasonable  assumptions  about 
A  dependence  on  particle  energy.  Thus  even  exponential  decays,  which  are  most 
probably  related  to  convection  and  adiabatic  deceleration  rather  than  to  diffu¬ 
sion,  are  in  some  cases  evidently  are  influenced  by  diffusion.  The  most  surprising 
result  is  the  presence  of  the  group  with  negative  a.  This  group  includes  events 
with  practically  all  possible  values  of  r  from  5  to  30  hours.  Since  the  background 
was  subtracted,  this  means  that  a  negative  (a  is  not  simply  an  instrumental  ef¬ 
fect  and  might  indicate  that  some  additional  exotic  effects  of  propagation  and/or 
acceleration  can  play  a  role. 

4.  Dependence  of  characteristic  decay  time  on  the  flare  site 

If  the  parent  flare  occurred  near  WOO  heliolongitude,  then  particles  accel¬ 
erated  at  the  flare  site  would  arrive  to  the  observer  along  the  shortest  way.  If 
the  flare  took  place  at  other  solar  meridian  one  could  expect  another  rise  time 
and  slower  decay  after  peak  intensity.  The  properties  of  decay  phase  of  a  few 
solar  flare  events  according  to  Pioneer  6-9  observations  was  studied  in  [5].  They 
obtained  that  in  events  registered  simaltaneously  at  different  heliolongitudes  r 
differs  by  l/^'o  d‘tp/dt  where  d'il)/dt  =  0.54  deg  hr tpo  is  about  ±30  deg.  Indi¬ 
rect  indication  for  such  a  dependence  was  also  obtained  in  [3],  where  AT  was 
considered  as  a  function  of  V’-  We  examined  if  such  dependence  exists  for  our 
whole  statistics  of  events.  Based  on  more  than  200  events,  the  dependence  r  on 
the  parent  flare  heliolongitude  during  1974-2001  was  analyzed.  The  result  is  that 
r  is  statistically  independent  from  the  flare  heliolongitude  (Figure  3),  while  the 
particle  arrival  and  peak  time  in  some  cases  are  significantly  delayed  with  the 
increase  of  the  distance  of  the  flare  site  from  the  optimal  heliolongitude  where 
magnetic  field  line  is  connected  to  the  observation  point. 

In  [5]  indications  were  also  obtained  for  the  events  considered  that  at  late 
times  (>  2  days)  the  spectral  exponent,  7,  near  10  MeV  depends  on  the  heliolon¬ 
gitude  of  the  observer  relative  to  the  centroid  of  the  particle  population  injected 
by  the  flare.  This  effect  results  in  a  variation  in  spectral  exponent  over  the  range 
2.0  <  7  <  4.5.  Again,  our  statistical  consideration  does  not  confirm  this  ob¬ 
servation.  This  suggests  that  after  SEP  event  maximum  particle  propagation  in 
the  interplanetary  medium  is  controlled  mainly  by  IM  conditions  in  the  neighbor¬ 
hood  of  observation  point,  and  the  structure  of  coronal  magnetic  fields  does  not 
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Fig.  3.  Dependence  of  r  and  a  on  heliolongitude  for  4.6-15  MeV  protons. 


influence  on  the  shape  of  time  profile  during  this  phase. 

5.  Conclusions 

1)  The  distributions  of  total  durations,  AT,  of  exponential  declines  with  and 
without  shocks  diflfer  insignificantly  for  AT  >  10  hrs,  while  for  AT  <  10  hrs 
the  relative  number  of  shock  associated  events  is  30-40  %  higher. 

2)  The  distribution  of  exponent  cx  firom  r  =  CE~°‘  representation  consists  of 
3  separate  groups:  i)  r  independent  of  proton  energy  (—0.1  <  a  <  0.1,  40 
events);  ii)  r  decreasing  with  energy  (a  >  0.1,  57  events),  and  iii)  a  small 
group  of  T  increasing  with  energy  (a  <  -0.1,  12  events). 

3)  Both  the  characteristic  decay  time,  r,  and  spectral  index,  7,  are  statistically 
independent  of  the  parent  flare  heliolongitude. 
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